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I review recent results derived from numerical simulations of the turbulent inter- 
stellar medium (ISM), in particular concerning the nature and formation of turbu- 
lent clouds, methods for comparing the structure in simulations and observations, 
and the effects of projection of three-dimensional structures onto two dimensions. 
Clouds formed as turbulent density fluctuations are probably not confined by ther- 
mal pressure, but rather their morphology may be determined by the large-scale 
velocity field. Also, they may have shorter lifetimes than normally believed, as 
the large-scale turbulent modes have larger associated velocities than the clouds' 
internal velocity dispersions. Structural characterization algorithms have started 
to distinguish the best fitting simulations to a particular observation, and have 
opened several new questions, such as the nature of the observed line width-size 
relation and of the relation between the structures seen in channel maps and the 
true spatial distribution of the density and velocity fields. The velocity field ap- 
parently dominates the morphology seen in intensity channel maps, at least in 
cases when the density field exhibits power spectra steep enough. Furthermore, 
the selection of scattered fluid parcels along the line of sight (LOS) by their LOS- 
velocity inherent to the construction of spectroscopic data may introduce spurious 
small-scale structure in high spectral resolution channel maps. 



1 Introduction 

Turbulence is a prime example of a chaotic system, and the interstellar 
medium (ISM) is most probably a prime example of a turbulent medium.E3 A 
discussion of interstellar turbulence thus befits this volume. Although chaos 
theory generally refers to systems with only a few degrees of freedom while 
turbulent flows have in general an extremely large number of them, both 
types of systems share the properties of sensitivity to initial conditions and 
the resulting practical unpredictability, as a consequence of the nonlinear cou- 
plings between the relevant variables. Furthermore, interstellar turbulence is 
much more complex than natural terrestrial and laboratory turbulence be- 
cause the former is magnetized, and, in cloudy regions, highly cortmressible 
and strongly self-gravitating, thus not expected to Kolmogorov-like,L3 except 
possibly in the diffuse gas. 

In recent years, many reviews covering various aspects of interstellar and 
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molecular cloud turbulence have appeared in the literature. Compressible tuL^ 
bulence basics and self-similar models are discussed by Vazquez-SemadeniH3 
A compendium of a wide variety of interstellar turbulence aspects is given 
in the volume Interstellar TurbulencejEi including in .particular turbulence in 
the HI gagj and in the diffuse ionized component .EjEjS A thorough review of 
the implications of compressible MHD turbulence fcaimolecular cloud and star 
formation is presented in Vazquez-Semadeni et al.Ej The present paper may 
be regarded as a companion to the latter reference, as it includes a number of 
topics not covered there. After reviewing the scenario of interstellar clouds as 
turbulent density fluctuations and some of its implications (sec. |2|), I discuss 
the comparison of the structural properties of simulated clouds with those 
derived observationally for real clouds (sec. H see also Ossenkopf et al., this 
volume, for a discussion from a more observationally-oriented perspective), 
and on the effects of projecting three-dimensional (3D) data onto two dimen- 
sions, the latter point being important for the interpretation of observational 
data, which necessarily are projections on the plane of the sky (sec. 

2 Interstellar Clouds as Turbulent Density Fluctuations 

A fundamental problem in the understanding of star formation is how the gas 
transits from a low-density diffuse medium to a comparatively enormously 
denser star. An intermediate step in this process is the formation of what we 
may generically refer to as interstellar clouds, including structures that span 
a wide range of physical conditions, from large diffuse HI clouds of densities a 
few X 10 cm"'^ and sizes up to hundreds of parsecs, to molecular cloud cores 
with densities ^ 10* cm"-^ and sizes of a few x 0.01 pc. 

The process of cloud formation quite possibly involves more than a sin- 
gle mechanism, including the passaEe.pf spiral density waves and the effects 
of combined large-scale instabilitiesEjO operating preferentially in the for- 
mation of the largest high-density structures, and-.the production of smaller 
density condensatipiis_.bjf jjiither swept-up shells Ji3'EllE3E3 or by a generally 
turbulent medium. BaoO'E^ In the remainder of this section we focus on the 
latter process. 

An important question is whether structures formed by either turbulent 
compressions or passages of single shock waves can become gravitationally 
unstable and collapse. This depends crucially on the cooling ability of the 
flow, which, as a first approximation, can be parameterized by an effective 
polytropic exponent jcS such that the pressure P behaves as P cx p'^"" , where 
p is the fluid density. Note that in this description both the "cooling" and 
the "pressure" can be generalized to refer to non-thermal energy forms, such 
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as magnetic and turbulent. 

The production and statistics of the density fluctuatkiLji&.in polytropic 
turbulence has beepJnvestigated recently by various groups.L§E3E3 Passot and 
Vazquez-Semadenilj have found that the probability density function (PDF) 
of the density fluctuations depends differently on the Mach number and on 
7efi. By means of a simple phenomenological model, these authors find that, 
for isothermal flows {'jcS = 1), the PDF is lognormal, as a consequence of the 
Central Limit Theorem and of the accumulative and multiplicative (additive 
in the log) nature of the density jumps caused by shocks. Increasing the 
rms Mach number only increases the width of the lognormal PDF and shifts 
its peak towards smaller densities. Instead, varying 7off changes the form 
of the distribution, which develops a power law tail at large densities for 
< 7cff < 1, and at small densities for jcS > 1- This effect is a consequence 
of the modification of the lognormal by the local variation of the sound speed, 
which in the general polytropic case varies with the density as p(Trff~i)/2. 
Essentially, in the power-law side of the PDF the density fluctuations are 
dominated by the nonlinear advection term in the momentum equation, with 
an increasingly negligible contribution from the pressure at increasingly large 
(7off < 1) or small {-fcs > 1) densities, due to the decreasing sound speed, 
while on the opposite side of the PDF the pressure dominates, impeding large 
excursions of the density. A concise disciiasion of the mechanism can also be 
found in Vazquez-Semadeni and Passot.E^ 

The stability of fluid parcels compressed, in n dimensions by shocks or 
turbulence requires 7off > 7cr = 2(1 — l/n)Hj'E3 For three-dimensional com- 
pressions, the minimum Mach number Mq necessary to induce collapse by 
the velocity field has been computed by several authors as a function of 7off 
and the mass m of the cloud in units of the —Leans mass. It is found that 
Mo oc In TO forJJie isothermal (7eff = 1) case ,Ea Mo oc TO('T=«'zi)/(4-37ctt) fo,- 
4/3 > 7eff > 1(N) and Mo > VlO/3(l-7eff) for < 7off < iM 

A relevant implication is that clouds formed by turbulent compressions 
are by necessity of a dynamical character, and are expected to either col- 
lapse, if the conditions described in the previous paragraph are satisfied, or 
else should reach some maximum density and then "rebound" as the external 
compression subsides (since the turbulent motions are chaotic, a compression 
will in general last a finite time only). This result has interesting implications 
on two "canonical" conceuts of interstellar dynamics, as recently discussed by 
Ballesteros-Paredes et al.u First, it may be that clouds JiL±he ISM may not 
be pressure-confined as in popular models of the ISM,l3E3 but rather in a 
highly dynamic and transient state, except if they become strongly gravita- 
tionally bound. An interesting corollary of this scenario occurs for regimes in 
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which the thermal jcS ^ (i.e., a nearly isobaric behavin^j, as is indeed the 
case for the ISM between densities '-^0.1 and 10^ cni~^.Ell In this case, the 
thermal pressure remains nearly constant, regardless of the density structures 
formed by the turbulence. However, while traditionally this near pressure 
constancy has been regarded as a pressure balance condition that provides 
confinement for clouds, in the turbulent-cloud scenario it is only a relatively 
irrelevant consequence of the medium being maintained at nearly constant 
thermal pressure by the prevailing cooling processes as it is compressed by 
the turbulent motions. 

Secondly, it appears difficult to produce the quasi-hydrostatic clumps 
which are the corninonly assumed to be the initial conditions of many models 
of star formation.Lj Their formation within a globally gravitationally stable re- 
gion by a turbulent compression requires, as described above, that 7off < 7cr- 
However, this sets them in a state of gravitational collapse. To then form 
a hydrostatic structure, a change in 7eff is required during the collapse, in 
order for to now become > 7cr, so that the pressure may now overcome 
the ongoing gravitational compression, and ultimately halt the collapse. Such 
change in 7cff, at least from thermal contributions alone, is not expected until 
very high densities ( ^ 10^ cm^^) are reached. □ 

Another implication of the turbulent density fluctuation scenario for the 
clouds is that the time scales associated to clouds may be smaller than those 
derived from their sizes and internal velocity dispersions.El If the cloud is 
made by the collision of large-scale streams, the time scale for its formation 
is the crossing time through its size L at the velocity difference between the 
turbulent scales larger than cloud (the colliding streams). For turbulence 
with a (normal) spectrum that decays with wavenumber fc, the characteristic 
velocity difference increases with separation, and thus the crossing time scale 
for the cloud is smaller than that derived from its internal velocity dispersion. 
Indeed, the line widths from the HI envelopes of molecular cloudS|are generally 
larger than the line widths in the molecular clouds themselves.O This result 
has been proposed as a possible solution to the absence of post-T-Tauri stars 
in the Taurus region, since the region may be younger than the age derived 
from its internal velocity dispersion.Q 



3 Comparisons Between Simulations and Observations 

In recent years, numerical simulations of interstellar turbulence have advanced 
to the point that statistical comparisons with observational results have be- 
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come possible.|j This is a crucial task because it allows an iterative procedure 
in which simulations may be constrained as models of the ISM and interstel- 
lar clouds by comparison of their morphological, topological and statistical 
properties with their observational counterparts. Once the best-fitting set of 
parameters is found for a certain type of system, the simulations may then be 
used as highly complete models of such system to improve our understand- 
ing of the physical processes occurring within them. However, it should be 
pointed out that this is a difficult task, because the natures of the observations 
and of the simulations are quite different. While simulations are performed 
on a regular grid, with well-defined boundaries, observations refer to regions 
of space for which the size along the line of sight is not constant. For example, 
for HI observations, the path length through the disk decreases with Galactic 
latitude, while for molecular line observations, the observed objects in general 
may have different extents along the LOS over their projected area. Also, 
linear sizes perpendicular to the LOS increase with distance. In spite of these 
difhcultics, however, several first steps have been taken in this direction. 

3.1 Scaling Relations 

One basic question-is whether the clouds in the simulations reproduce the 
well-known LarsonEa relations (see also the review by BlitzQ for a more recent 
account) Av oc i?^/^ and (p) oc , where Av is the velocity dispersion in the 
cloud, (p) its mean density, and R its characteristic size. In an analysis of two- 
dimensional MHD simulations of the ISM including self-gravity and stellar-like 
driving, Vazquez-Semadeni et al.E3 found that, although with much scatter, 
a Larson-like velocity dispersion-to-size scaling of the form At; cx is 
observed for clouds defined as the connected regions in the flow with densities 
above a given threshold. This result is roughly consistent with observational 
surveys giving scaling exponents between 0.4 and 0.7.0 However, the density- 
size relation is not verified in the simulations. Instead, small clouds with 
low densities, which are transient and not gravitationally bound, are formed 
in large quantities in the simulations. Rather than being satisfied by all 
clouds, the Larson density-size relation appears to be an upper bound to 
the region populated by the clouds in a p-R diagram. The same treiud was 



the limitations on integration times of observational surveys, and that the 
Av-R relation, satisfied by all clouds, may originate from the Burgers-like 

"^But note that modeling of individual objects is not feasible because of the sensitivity to 
initial conditions of turbulent flows. 
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spectrum of molecular cloud turbulenceO In this scenario, only those clouds 
which become sufficiently self gravitating and are not strongly disturbed (i.e., 
in near virial equilibrium), satisfy the density-size relation. 

3.2 Synthetic Line Profiles 

Another important means of comparison between simulations and observa- 
tions are the spectral line profiles, since their shapes (generally character- 
ized by their first few moments) reflect the velocity structure of the flow in 
the observed regions. Line profiles from the simulations are constructed as 
density-weighted vebcity histograms along each LOS. 

Falgaronejet alEj compared the line spectra produced in a high-resolution 
3D simulationEll of weakly compressible turbulence with observational data, 
concluding that both sets of spectral lines are very similar, in terms of the 
range of values of the variance and the kurtosis they present. The similarity 
is greater in this case than that achieved by other models of clouds (con- 
structed with random uncorrelated velocity fields or with isolated clumps in 
an interclump medium) which do not account for the spatial correlations de- 
rived from the continuum nature of fiuid turbulence. Similar results have 
been obtained from randomly generated flows with an imposed Kolmogorov 
spectrum.EJ However, a more recent studx based on the PDFs of the velocity 
centroids, rather than on the line proflleSjCHI suggests that in turn nearly incom- 
pressible turbulence fails to capture some features of the centroid PDFs of both 
molecular and HI regions, which exhibit a larger degree of non- Gaussian! ty 
than those derived from incompressible or weakly compressible turbulence. 

Line profiles from simulations of strongly compressible MHD turbulence, 
including radiative transfer, as well as other diagnostics, have recently been 
compared to molecular line data by Padoan and coworkers to support their 
suggestion that molecular clouds may actually be in a. super- Alfvenic regime, 
rather than a sub- Alfvenic one, as generally believedfl Their tests were based 
on two 3D simulations of MHD isothermal turbulence without self-gravity, one 
super- Alfvenic, the other sub- Alfvenic. First, they noted that simulated line 
profiles from the super- Alfvenic run seem to reproduce the observed growth of 
line width with integrated antenna temperaturec^ better than those from the 
sub- Alfvenic run.O Secondly, the super-Alfvenic run seems toJpettcr match 
the observed trend of magnetic field strength vs. gas density.^ Finally, the 
super-Alfvenic simulation also more closely reproduces the observational trend 
of the_jdispersion of extinction vs. mean extinction along selected lines of 
sight.cS However, one caveat remains before the super-Alfvenic model can be 
accepted: the simulations lacked self-gravity, which could have pushed the 
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results of the sub-Alfvenic simulation closer to the observational results, and 
the super-Alfvenic run away from them. Therefore, similar experiments are 
needed with self-gravitating runs in order to confirm this possibility. 



3.3 Higher-Order Statistics and Fractality Analyses 

The methods described in the previous section have taken into account only 
the velocity information in the spectral data "cubes" . We now briefly discuss 
methods that have tried to take spatial information into consideration as well. 

Spatial structure is often described by means of the autocorrelation func- 
tion (ACF), which measures the probability of finding equal values of a given 
physical variable at two different positions in space, as a function of their sep- 
aration.^ Early studies in this direction measured the ACF of column density 
and of the line velocity centroid&^ajttenmting to determine whether charac- 



teristic lengths exist in the IS 



with mixed results. Recently, a 



variant of this approach termed the Spectral Correlation Function (SCF) has 
been introducedJla The SCF measures the quadratic difference between line 
spectra (on a channel-by-channel basis) at different positions on a spectral- 
line map, in an attempt to include spectral as well as spatial information in 
the statistical description. So far, the method has been used to measure the 
angle-averaged correlation between the spectrum at a given position in a map 
and at its nearest neighbors, allowing the characterization of the small-scale 
variability of the spectra, and a comparison between CO maps and simula- 
tions of isothermal turbulence under various regimes (purely hydrodynamic, 
MHD, and self-gravitating). In that work, differences between theo^alues of 
the SCF for weakly compressible purely hydrodynamic turbulenceEll and for 
the Ursa Major molecular cloud indicate that simulation is not as accurate_a 
model for the Ursa Major cloud as previously claimed by Falgarone et alc3 
on the basis of line profile shapes (see sec. |3.2| ). Comparison of HI data with 
non-isothermal simulations has proven more difficult, because the thernial 
broadening of the warm gas swamps the velocity structure of the cold gas.Q 

The recognition that the structure of the ISM may have a turbulent ori- 
gin has also prompted searches for fractal properties of interstellar clouds 
(although see the contribution by Combes, this volume, for an alterna±ii£fi 
scenario originating the fractal structure). Early studies in this directionEilO 
started by measuring the fractal dimension of the clouds in column density 
or intensity maps of selected regions, by means of the area-perimeter scaling 

''However, it has been pointed out by Scale and Houlahan l^F^ that one limitation of the 
ACF is that it cannot distinguish between hierarchically nested or randomly distributed 
structure. 
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in the projected clouds, finding dimensions near 1.4. Recent measurements of 
their area-perimeter relation in 2D simulations find similar values. Q However, 



these values are surprising 
in the Earth's atmosphere 



close to the projected fractal dimensions of clouds 
I, suggesting that the fractal dimension may be a 
significantly degenerate diagnostic which may not be capable of distinguishing 
between different physical regimes. Indeed, a more recent study by Chappell 
and Scalollj has shown that the column density maps of various regions actu- 
ally have a well-defined multifractal structure, so that attempting to measure 
a single fractal dimension for the clouds in the maps may erase much of the 
structural information. Furthermore, these authors emphasize that the meth- 
ods used to determine fractal dimensions of clouds rely on the definition of 
"clouds" by means of some rather arbitrary criterion (such as thresholding 
the column density field) , while the multifractal spectrum determination uses 
the structural information of the whole field. The multifractal spectrum of 
the regions studied appears to correlate fairly well with the geometric forms 
seen visually, potentially providing a means for quantitative structure classi- 
fication schemes. The multifractal propertiejS_pf numerical simulations of ISM 
turbulence are currently being investigated.Ej 

A method for determining the line width-size scaling of spectral maps 
independently of any specific definition of "clouds" pii-ai spectral map has 
been introduced recently by Heyer and coUaboratorsBJEa The method uses 
the statistical technique known as Principal Component Analysis (PCA) to 
define a set of spectral profiles (eigenvectors) which form a "natural" basis for 
the spectral maps (in the sense that it reflects the main trends of the intensity 
data among the velocity channels). Eigenimages, which represent the intensity 
structures as "filtered" by the basis spectra, are generated by projecting the 
original spectra onto the eigenvectors. By then measuring the ACF of the 
eigenvectors and of the eigenimages, the relationship between the magnitude 
of velocity differences and the spatial scales over which these differences occur 
can be extracted. In order to "calibrate" the expected value of the exponent 
a in this relation (such that At; cx i?", where Aw is the velocity difference 
and R is the size) , "pseudo-simulations" of fractional Brownian noise with a 
prescribed-Spectrum were produced and "observed" using a radiative transfer 
simulator .c3 It was found that the scaling exponent is related to the spectral 
index_/3 by a = 13/3. Recent calibrations with actual hydrodynamic and 
MHDc3 simulations in 3D are roughly consistent with this result, although 
the resolutions available in 3D simulations are still insufficient to develop 
clear power-law turbulent spectra, so the results are not conclusive. Tests with 
higher resolution 2D simulations apparently produce a different calibration, 
a = (3/4:. 
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The origin of these scahngs is not well understood yet. In fact, they 
suggest that the very nature of the line width-size relation derived through 
PCA is unknown, since it does not follow the same scaling as, say, the second- 
order structure function, defined as F2{r) = ([u(x) — u(x + r)]^), where the 
brackets denote a volume average. This function gives the mean quadratic 
velocity difference betweenjiositions separated a distance r. Yet, F2 is related 



so that, if E{k) oc k~^, then F2{r) cx , with rj = {(3 — l)/2. 7] is thus not 
functionally related to (3 in the same manner as a. The discrepancy is probably 
related to the fact that the above scaling refers to the structure function of 
the actual 3D velocity field, while in spectroscopic data every velocity interval 
contains the contribution of many fluid parcels (possibly far apart from each 
other), and only one of the three velocity components is observed. Thus, 
the true nature of the spatial and velocity increments in spectral data cubes 
remains unknown. 

Finally, a method for structure analysis similar to the power spectrum, 
but using a wavelet rather than a Fourier basis is described in the contribution 
by Ossenkopf et al. (this volume). 

4 Effects of Projection on Morphology 

One advantage of 3D simulations is that they contain more structural infor- 
mation than even spectroscopic data "cubes" . While the latter only refer to 
two spatial and one velocity dimensions, 3D numerical simulations provide 
information on the 3D structure of all physical variables. (Of course, their 
downside is that they are necessarily limited in resolution and in the number 
of physical processes that can be included.) This allows an investigation of 
the 3D structures that generate the patterns observed in the position-position- 
velocity (PPV) space of the spectroscopic channel maps. To this end, channel 
maps are constructed from the simulations by integrating the density fleld 
along one direction (the line of sight, or LOS), and then selecting the contri- 
bution of each parcel along the LOS by its LOS-velocity. This is equivalent 
to constructing density- weighted velocity histograms (the "line spectra") at 
each position in the plane perpendicular to the LOS. 

A rather unexpected result has recently been found independently, using 
different approaches, by several groups. It appears that, at least under certain 
conditions, the projected spatial structure in the channel maps is dominated 
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to the energy spectrum 




(1) 



by the spatial distribution of the velocity field, rather than by the 3D density 
field. Pichardo et al.E3 have shown this in a 3D simulation of the ISM at inter- 
mediate scales (3-300 pc) by noting that the pixel-to-pixel correlation between 
channel maps and thin slices of the 3D velocity field tends to be larger on av- 
erage than the correlation between channel roaps and slices of the 3D density 
field. Independently, Lazarian and PogosyantJ have shown analytically that, 
for cases with an underlying one-to-one mapping between the position along 
the LOS and the LOS-velocity (as for an expanding universe or the HI gas dis- 
tribution in the Galaxy), and with uncorrelated random density and velocity 
fields with well-defined spectral indices, the power spectrum of the projected 
density field is dominated by the spectrum of the velocity field for density 
spectral indices steeper than —3, unless the velocity channels are very wide 
(as is clearly the case in the limit of a single velocity channel, m_which the 
velocity dependence is integrated out). Finally, He yer and BruntSS have no- 
ticed that, in their pseudo-simulations (cf. sec. |3.3| above), channel maps and 
PCA-derived Av-R relations produced with and without density-weighting 
are notoriously similar, suggesting that the effect of the density weifihting is 
relatively minor. A similar effect was noticed by Falgarone et al.EJ about 
the shapes of synthetic line profiles. In summary, it appears that the spatial 
structure of the velocity field is at least as important as that of the density 
field in determining what is observed in projection on theufllane of the sky. 

A related effect has been observed by Pichardo et al.cS The morphology 
observed in the channel maps appears to contain much more small-scale struc- 
ture than either the density or the velocity 3D fields. This is reflected in the 
power spectra of the channel maps and of 2D slices through the 3D density and 
velocity fields, the latter two having steeper slopes and falling off much more 
rapidly than the former. This phenomenon has been interpreted by those au- 
thors as a consequence of a pseudo-random sampling of fluid parcels along the 
LOS by the velocity selection performed when constructing a channel map. 
This introduces an additional ingredient of variability between neighboring 
LOSs, which causes artificial small-scale variability in the channel maps. 

It can be concluded from this section that, for a fully turbulent ISM, the 
structure seen observationally, through spectroscopic observations, may differ 
from the actual 3D structures present in the medium. In particular, this 
suggests that structure-finding algorithms operating on spectroscopic data 
cubes may not identify exactly the same structures than would be obtained on 
the actual 3D spatial data, as already pointed out by various authors. Ella These 
effects may be decreased, however, in cases when the observed regions contain 
well-defined "objects" which may be picked out by the observing process, such 
as shells, bipolar flows, etc. 
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5 Conclusions 



In this review I have discussed recent resuhs from numerical simulations of 
turbulence in the ISM. I first reviewed the scenario of interstellar clouds as tur- 
bulent density fluctuations. Work on the production of gravitationally bound 
structures in globally stable media by turbuleiit, compressions was summa- 
rized, in particular the necessary Mach numberso (for 3D compressions) and 
the constraints, on the effective polytropic exponent 7off|-,for n-dimensional 
compressions. e3 Three implications were then discussed u'EI First, the near 
thermal pressure balance observed in the ISM (except for molecular clouds) 
may not be a confining agent for clouds, but rather a relatively fortuitous 
consequence of the prevailing heating and cooling mechanisms, which render 
the medium nearly isobaric, in the presence of turbulence-induced density 
fluctuations. Second, it appears unlikely that nearly hydrostatic cores may be 
produced within the turbulent ISM unless some very specific variations in 7efT 
occur during a gravitational contraction induced by the turbulence. Third, 
the time scales associated with clouds may be smaller than those estimated 
from the clouds' characteristic dimensions and their velocity dispersions, since 
the relevant velocities may instead be those of the larger, external flow streams 
that produced the clouds at their collision interfaces. 

I then proceeded to review recent results from various attempts to relate 
numerical simulations to ohaervational data, from early qualitative compar- 
isons of spectral line profilea£d and surveys of clouds in 2D simulations (which 
suggested the existence of a whole population oOow-column density clouds 
that do not satisfy Larson's density-size relation)E3, to recent approaches us- 
ing morcj^phisticated statistical techniques s.uch as the Spectral Correlation 
Functioned, Principal lOpmponent AnalysiscSE3, and fractal dimensions and 
multifractal spectraE2l'E2lll3, mostly aiming at characterizing the morphology 
of interstellar structures in a statistically meaningful way, and determining 
whether the structures developing in turbulence simulations reproduce their 
properties. Some of these methods are only being developed now, but they 
are already providing a quantitative method for discriminating between turbu- 
lence simulations with different parameter choices as the most suitable models 
for specific interstellar regimes, as well as providing a basis for interpreting 
observational data in terms of the simulations. 

Finally, I discussed the relationship between the actual 3D spatial struc- 
ture of the density and velocity fields, and that of the projected-liiiitensity'' 
field in channel maps. Recent works, using various approachescJ'c£rE3, suggest 
that the structure, both morphological and statistical (power spectrum) of 
the 2D intensity field is dominated by the velocity field, rather than by the 
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density. Additionally, since the total intensity in every LOS of a channel map 
is constructed by "selecting" scattered fluid parcels (on the basis of their LOS 
velocity) along the LOS, spurious small-scale variability is introduced into 
the structure seen in the channel maps, as the set of sampled parcels varies 
from one LOS to the next.E^ These results suggest that the structure in chan- 
nel maps bears a very complex and non-trivial relationship to the structures 
actually existing in the ISM. Further work in this area is likely to produce 
numerous unexpected and exciting results in the near future. 
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